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Theoretical  analyses  are  performed,  and  numerical  results  obtained  for  the 
axisymmetric  scattering  of  a spherical  wave  by  either  a rigid  or  soft 
(pressure  release)  prolate  spheroid.  In  the  numerical  calculations,  the  com- 
plex far-field  reflection  factor  (ratio  of  scattered-to-incident  pressure) 
is  computed  as  a function  of  frequency  and  observation  angle  for  a 7-to-l 
fineness  ratio  prolate  spheroid.  Ten  different  source  locations  are  con- 
sidered, all  of  which  lie  on  the  axis  of  symmetry.  In  addition,  the  super-  - 
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20.  position  of  the  incident  and  scattered  fields  is  investigated  for  a 
variety  of  two-source  configurations.  In  one  instance,  these  two 
sources  are  assumed  to  have  equal  strengths  and  to  vibrate  in  phase, 
while  in  a second  instance  they  are  assumed  to  be  incoherent. _ 


The  computations  have  been  performed  for  relatively  low  frequencies, 
where  the  exact  tabulated  forms  of  the  spheroidal  wave  functions  may  be 
used.  A reduced  frequency,  defined  by  kd/2  ranged  from  0.1  to  40,  where 
k is  the  wavenumber,  and  d is  the  interfocal  distance  of  the  7-to-l 
spheroid.  The  results  for  the  directivity  functions  indicate  that  the 
reflection  factor  depends  only  slightly  on  aspect  angle  for  a single 
source,  but  when  two  in-phase  sources  are  considered,  the  scattering  is 
very  dependent  on  angle.  If  the  two  sources  are  incoherent,  then  the 
patterns  approach  the  shape  of  the  corresponding  pattern  for  the  single 
source  closest  to  the  pole,  particularly  when  the  second  source  is 
relatively  far  from  the  first.  The  reflection  factor  is  shown  to  increas 
with  frequency  in  the  rigid  case  and  to  remain  essentially  constant  with 
frequency  in  the  soft  case. 
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ABSTRACT 

Theoretical  analyses  are  performed,  and  numerical  results  obtained 
for  the  axisymmetric  scattering  of  a spherical  wave  by  either  a rigid  or 
soft  (pressure  release)  prolate  spheroid.  In  the  numerical  calculations, 
the  complex  far-field  reflection  factor  (ratio  of  scattered-to- incident 
pressure)  is  computed  as  a function  of  frequency  and  observation  angle  for 
a 7-to-l  fineness  ratio  prolate  spheroid.  Ten  different  source  locations 
are  considered,  all  of  which  lie  on  the  axis  of  symmetry.  In  addition, 
the  superposition  of  the  incident  and  scattered  fields  is  investigated  for 
a variety  of  two-source  configurations.  In  one  instance,  these  two  sources 
are  assumed  to  have  equal  strengths  and  to  vibrate  in  phase,  while  in  a 
second  instance  they  are  assumed  to  be  incoherent. 

The  computations  have  been  performed  for  relatively  low  frequencies, 
where  the  exact  tabulated  forms  of  the  spheroidal  wave  functions  may  be 
used.  A reduced  frequency,  defined  by  kd/2  ranges  from  0.1  to  40,  where 
k is  the  wavenumber,  and  d is  the  interfocal  distance  of  the  7-to-l 
spheroid.  The  results  for  the  directivity  functions  indicate  that  the 
reflection  factor  depends  only  slightly  on  aspect  angle  for  a single 
source,  but  when  two  in-phase  sources  are  considered,  the  scattering  is 
very  dependent  on  angle.  If  the  two  sources  are  Incoherent,  then  the 
patterns  approach  the  shape  of  the  corresponding  pattern  for  the  single 
source  closest  to  the  pole,  particularly  when  the  second  source  is 
relatively  far  from  the  first.  The  reflection  factor  is  shown  to  increase 
with  frequency  in  the  rigid  case  and  to  remain  essentially  constant  with 
frequency  in  the  soft  case. 
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CHAPTER  I 


INTRODUCTION 


1.1  General  Considerations 

The  problems  related  to  spheroids  are  of  various  types  like  diffraction 
and  scattering,  perturbation  in  a steady  flow,  radiation,  and  cavitation. 
They  are  commonly  encountered  in  acoustics,  as  well  as  in  hydrodynamics  or 
electromagnetics.  The  very  wide  variety  of  shapes  that  a spheroid  may 
take,  from  a sphere  to  a needle,  in  the  case  of  prolate  spheroids,  and 
from  a disk  to  a sphere,  in  the  case  of  oblate  spheroids,  makes  their  use 
very  practical. 

In  the  field  of  acoustics,  we  are  interested  in  radiation,  scattering, 
and  diffraction  due  to  spheroids.  For  that  purpose,  we  have  to  solve  the 
equation  for  wave  propagation.  As  the  spheroidal  coordinate  system  is  one 
of  the  eleven  coordinate  systems  in  which,  for  harmonic  time  dependence, 
the  reduced  wave  equation  (the  Helmholtz  equation)  is  separable,  the 
solution  to  the  wave  equation  can  be  expressed  in  terms  of  spheroidal 
eigenfunctions,  or  spheroidal  wave  functions. 

Mathematicians  and  physicists  have  investigated  the  spheroidal  wave 
equations,  in  order  to  find  expressions  leading  to  numerical  values  or 
asymptotic  expansions.  Between  1830  and  1860,  Lam6  first  developed 
solutions  to  the  Helmholtz  equation  in  general  ellipsoidal  coordinates. 
However,  their  complexity  did  not  allow  efficient  use.  Niven  (1880), 
noting  that  the  spheroidal  coordinate  systems  reduced  to  the  spherical 
coordinate  system  at  large  distances  from  the  origin,  expanded  the 
spheroidal  wave  functions  in  terms  of  the  spherical  ones.  Power  series 
expansions  were  introduced  by  MacLaurin  (1898).  Abraham  (1899)  obtained 
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some  integral  solutions;  those  were  treated  more  extensively  by  Poole 
(1923).  Then  Meixner  and  ShMfke  (1954)  performed  extensive  work  on 
spheroidal  wave  functions,  introducing  asymptotic  expansions.  Flammer 
(1957) made  a collection  of  the  preceding  works  and  provided  tables  for 
various  coefficients,  eigenvalues  and  eigenfunctions.  Silbiger  (1961) 
and  Slepian  (1965)  have  investigated  asymptotic  expansions  which  are  of 
great  use  today.  Recently,  Hanish  et  al.  (1970)  and  Van  Buren  et  al. 

(1975)  have  provided  extensive  tables  for  the  wave  functions  at  low 
frequencies. 

1.2  Problem  to  be  Investigated 

Different  efforts  have  taken  place  concerning  acoustic  radiation, 
scattering  and  diffraction  by  spheroids.  A very  extensive  bibliography 
concerning  spheroidal  wave  functions,  and  the  various  advancements  on  the 
acoustics  of  spheroids  can  be  found  in  Van  Buren  et  al.  (1975). 

In  this  study,  we  investigate  the  far-field  reflection  factor 
generated  by  a prolate  spheroid  of  fineness  ratio,  e = 7.1,  insonified  by 
a point  source  located  on  its  major  axis  at  a short  distance  from  the  >i 

pole.  Both  cases  of  a rigid  surface  and  of  a pressure  release  surface  are 

\ | } 

. 

[ 

i 


J 


investigated.  The  source  assumes  ten  different  positions.  Furthermore, 
we  consider  combinations  of  two  sources  in  phase  as  well  as  incoherent, 
one  source  being  fixed,  the  other  successively  occupying  the  remaining 
nine  positions.  For  all  these  configurations,  we  determine  the  directivity 
pattern  of  the  reflection  factor  for  the  values  of  the  reduced  frequency 
parameter:  h = kd/2  = 20,  35,  40,  where  k is  the  wavenumber,  and  d is 
the  interfocal  distance. 


We  determine  the  frequency  response  in  the  range  of  h from  0.1  to  20 
for  the  three  values  of  the  aspect  angle:  0=0°,  90°,  180°  as  measured 
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CHAPTER  II 


PROLATE  SPHEROIDAL  COORDINATES,  WAVE  FUNCTIONS, 
AND  GREEN'S  FUNCTION 


2.1  Introduction 

In  this  chapter  the  prolate  spheroidal  coordinate  system  is  discussed. 
The  Helmholtz  equation  is  expressed  in  that  coordinate  system  and  is  reduced 
to  two  ordinary  differential  equations  after  separation  of  variables. 

These  equations  are  solved,  leading  to  the  prolate  radial  and  prolate 
angular  spheroidal  wave  functions. 

The  approach  followed  to  solve  the  differential  equations  was  first 
introduced  by  Niven  (1880)  and  has  been  extended  by  Flammer  (1957). 

However,  the  normalization  scheme  followed  is  that  one  first  introduced 
by  Meixner  and  Shdfke  (1954)  and  used  by  Kanish  et  al.  (1970). 

Further,  the  derivation  of  the  free-space  Green's  function,  which  is 
of  fundamental  use  in  field  problems,  is  derived  in  terms  of  the  prolate 

spheroidal  wave  functions  according  to  the  scheme  given  by  Flammer  (1957). 

\ 

2.2  The  Prolate  Spheroidal  Coordinates 

The  prolate  spheroidal  coordinate  system  is  derived  from  the  general 
ellipsoidal  coordinate  system  by  choosing  the  major  axis  to  be  an  axis  of 
symmetry.  That  is,  it  may  be  generated  by  the  rotation  of  the  two- 
dimensional  elliptic  coordinate  system  about  its  major  axis  (Figure  1). 

The  alternative  possibility  of  rotating  the  two-dimensional  elliptic 
coordinate  system  wtih  respect  to  its  minor  axis  generates  the  oblate 
spheroidal  coordinate  system. 

An  ellipse  is  the  locus  of  points  such  that  the  sum  of.  the  distances 
r^  and  r ^ from  two  fixed  points  (the  foci),  separated  by  a distance  d,  is 
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1 


k 


a constant,  a property  which  can  be  written  as 


rx  + r2  « £d  , 


(2.2.1) 


where  £ is  the  nondimensional  elliptic  coordinate.  Therefore,  after 
rotation  of  the  ellipse,  £ represents  a given  ellipsoid  and  by  analogy 
with  the  spherical  coordinate  system  is  called  the  spheroidal  radial 
coordinate.  For  a given  finite  value  of  d,  £ varies  from  1,  where  it 
represents  the  segment  connecting  the  two  foci,  to  infinity,  where  it 
represents  a sphere  of  infinite  radius,  the  center  of  which  is  located  at 
the  origin  of  the  coordinate  system.*  Some  characteristic  dimensions  of 
the  ellipsoid  represented  by  £ can  be  easily  derived:  the  length  of  the 
spheroid  is  given  by 

L = ri  + r2  » £d  , (2.2.2) 


the  minor  radius,  R,  is  obtained  by 


2 ^ ,d. 2 22 

R + (j)  = rl  + r2 


and 


noting  that  r^  = = £d/2,  we  get 


R = | 42  - 1 5 


therefore,  the  fineness  ratio,  e,  is  given  by 


2R 


k1  - 1 


(2.2.3) 


(2.2.4) 


The  radius  of  curvature  R^  of  the  spheroid  varies  with  the 
observation  point  on  its  surface  and  can  be  obtained  by  looking  at  an 
elliptic  section  in  the  x-z  plane  of  Figure  1 and  writing  the  general 
relation 

o _ (1  + x'2)3/2 


*If  we  let  £ ■+  oo  while  d -*•  0,  keeping  £d  « D,  the  spheroid  tends  towards 
a sphere  of  finite  radius  D/2. 
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where  primes  denote  d/dz.  The  equation  of  the  ellipse  is  given  by 


2 2 
— + _1_ 

R2  (fo2 


Hence, 


16z2  L2  „2. . 3/2 
,4  l4  ' R )J 


(2.2.5) 


In  the  mid-section  the  extreme  values  of  the  radius  of  curvature 
are  R^  at  the  equator  and  R2  at  the  pole,  and  are  obtained  respectively 
by  letting  z = 0 and  z = L/2  in  Equation  (2.2.5);  we  obtain 

,2 


and 


R _ J±_ 

R1  “ 4R 


R -2R- 
R2  L 


(2.2.6) 


(2.2.7) 


The  spheroidal  coordinate  system  being  orthogonal  has  a second  family  of 
surfaces.  It  is  the  double  set  of  conjugate  hyperboloids  obtained  from 
the  same  focal  points  as  the  ellipsoids.  They  can  be  generated  by  rotating 
about  the  z axis  the  nonfocal  hyperbolas,  orthogonal  to  the  previous 
confocal  ellipses  (Figure  1).  Such  a hyperbola  is  the  locus  of  points 
such  that  the  difference  of  the  distances  r^  and  from  the  foci  is 


constant,  a property  which  can  be  written  as 


rl  ' r2  = nd 


(2.2.8) 


It  can  be  easily  seen  that  this  hyperbola,  for  large  values  of  r^  and 
is  asymptotic  to  a line  making  an  angle  cos  with  the  z axis.  That 
means  that  for  large  values  of  £ this  angle  is  identical  to  the  spherical 
aspect  angle,  i.e.. 


0 * cos 


(2.2.9) 
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Therefore,  n is  called  the  spheroidal  angle  coordinate.  As  each  value  of 
n defines  one  hyperboloid  and  only  one,  each  value  of  0 in  the  range  from 
zero  to  n also  defines  a unique  hyperboloid  asymptotic  to  the  cone  of  angle 
0 with  respect  to  the  z axis.  Therefore,  n must  range  from  1 to  -1  to 
cover  the  whole  family  of  hyperboloids.  When  n = 0,  it  refers  to  the 
equatorial  plane  of  the  spheroid,  the  hyperbola  degenerating  into  the 
x-y  plane.  When  n = +1 , it  refers  to  the  pole  region  located  on  the 
positive  side  of  the  z axis  and,  when  n ■ -1,  it  refers  to  the  pole  region 
located  on  the  negative  side  of  the  z-axis,  the  angle  0 being  here  defined 
from  the  positive  side  of  the  z-axis  (Figure  1) . 

The  third  family  of  surfaces  used  in  the  spheroidal  coordinate  system 
are  the  surfaces  orthogonal  to  both  the  ellipsoids  and  the  hyperboloids, 
i.e.,  the  planes  containing  the  z-axis.  As  the  previous  surfaces  have 
been  generated  by  rotation  of  plane  curves  about  this  axis,  the  third 
coordinate  $ represents  the  rotational  angle  and  is  called  the  rotational 
coordinate  (Figure  1).  It  ranges  from  zero  to  2ir  and  is  similar  to  the 
rotational  coordinate  in  both  spherical  and  cylindrical  systems. 

From  the  previous  remarks,  we  can  see  that  for  large  values  of  £, 
where  the  spheroid  reduces  to  a sphere,  the  spheroidal  coordinate  system 
reduces  to  the  spherical  coordinate  system. 

The  following  relations  between  the  Cartesian  coordinates  x,  y,  z 
and  the  spheroidal  coordinates  C,  n,  <P  can  be  easily  derived  and  are  given 
in  many  works  like  Flammer  (1957),  Hanish  et  al.  (1970),  Skudrzyk  (1971), 
and  Van  Buren  et  al.  (1975): 

x - |[(1  - n2)(C2  - 1 ) 1 1 /2  * cos  * , (2.2.10) 

y - f[(l  - n2)U2  - l)]1/2-sin  4>  , (2.2.11) 
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z = f nS 


(2.2.12) 


2.3  The  Wave  Equation  In  Prolate  Spheroidal  Coordinates 

As  our  Interest  is  in  investigating  wave  propagation  in  spheroidal 
coordinates,  we  first  solve  the  equation  of  propagation  in  that  system 
by  separation  of  variables. 


-Mf  - c2V2*  = 0 , 
at 


(2.3.1) 


where  t is  the  time  variable,  and  c the  wave  propagation  velocity.  As  we 
will  be  concerned  with  harmonic  time  dependence,  the  function  of  time  T 
reduces  to 


T = T e 
o 


(2.3.2) 


where  w is  the  angular  frequency. 

Our  investigation  will  then  be  concentrated  on  solving  the  remaining 
part  of  the  wave  equation,  representing  the  spatial  dependence  of  pressure 
(or  any  chosen  quantity  like  particle  velocity,  density , velocity  potential , 
temperature),  i.e.,  the  classic  Helmholtz  equation  is: 


(V2  + k2)i|/  = 0 , 


(2.3.3) 


where  ip  represents  a general  variable,  and 


k = c 


(2. 3. A) 


is  the  wave  number. 

The  Laplacian  operator  can  be  derived  by  using  the  formulas  for  the 
divergence  and  for  the  gradient  in  spheroidal  coordinates  given  in  Hanish 
et  al.  (1970),  and  Van  Buren  et  al.  (1975).  It  can  also  be  computed 


directly  from  the  metric  coefficients  as  in  Flammer  (1957)  and  Lauchle 
(1970).  This  leads  to: 
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V2  = 


2 2 

cr<r 


_fJ_„  2._8_  , 2 ._3_ 

2X  3r>  1 n + 3£(^  " 1J3£ 


n ) 


r2  2 
£ - n 


2 2 2 J 

(5  - DU  - n ) 3<T 


(2.3.5) 


After  substitution  into  Equation  (2.3.3)  and  using 

kd 


(2.3.6) 


as  a nondimens ional  frequency  parameter,  the  Helmholtz  equation  can  be 
written: 


r J_/i  2x_3_  . _3^  1V_3_  , 

[3n(1  “ n 3£(C  " X)3C  + -.2 


.2 

3 


<f,‘  - i)d  - n2)  3<t>2 


+ h2(C2 


n )]*  - 0 . (2.3.7) 
As  the  prolate  spheroidal  coordinate  system  is  orthogonal  we  can  use  the 
separation  of  variables  method  to  solve  that  equation.  Solutions  can  be 
obtained  in  the  form  of  Lamp's  products: 


L = smn(h*h)R  (li.O-T  H 

mn  mn  mn  sin 


(2.3.8) 


where  S (h,n)  and  R (h,£)  are,  respectively,  the  angle  and  the  radial 
mn  mn 


spheroidal  wave  functions.  Substituting  the  right-hand  side  of  Equation 
(2.3.8)  into  Equation  (2.3.7)  leads,  after  some  calculations  presented  in 
Lauchle  (1970),  to  two  ordinary  differential  equations: 


-r-[(l  - n2)-r-  S (h,n)]  + [X  - h2n2 
d p dp  mn  mn 


m 


]S  (h,n)  = 0 (2. 3. 8. a) 


1 - n 


21  u 

mn 


j 
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-£r[a2  - l)j=  R (h.O ] - U ~ h2C2 

□ 4 mn  mn 


- pTT"”<h'0 ' 0 • 


(2.3.8  .b) 


where  X are  eigenvalues  and  are  identical  for  both  equations. 

These  two  equations  appear  to  be  similar  and  are  of  the  type: 
a ~ du  2 

-^[(1  - z2)-^l  + [X  " hV — j]u  = 0 . (2.3.9) 

dz  dz  mn  , 2 mn 

1 - z 

The  solution  of  Equations  (2.3.8)  are  particular  ellipsoidal  wave 
functions,  or  Lamp's  wave  functions,  for  the  degenerate  case  of  the 
spheroidal  coordinates.  However,  they  are  difficult  to  determine  and 
not  useful  for  practical  applications.  Power  series  solutions  are  not 
of  great  use  either,  due  to  the  singularities  presented  by  Equation 
(2.3.9)  for  z = + 1 and  z ■+  »>.  Niven  (1880)  first  investigated  these 
equations  and  expressed  the  spheroidal  wave  functions  as  series  of  the 
spherical  wave  functions,  as  a consequence  of  the  identity  of  spheroidal 
and  spherical  coordinate  systems  for  large  values  of  £.  The  radial 
functions  are  expressed  in  terms  of  spherical  Bessel  functions  and  the 
angle  functions  in  terms  of  associated  Legendre  functions.  A detailed 
derivation  of  those  expressions  is  given  in  Flammer  (1957)  and  Lauchle 
(1970);  therefore,  we  simply  give  an  outline  here. 


2. A The  Angle  Functions 

We  have  to  solve  Equation  (2. 3. 8. a).  We  must  first  note  that  in  most 


physical  applications,  the  continuity  of  the  field  on  a line  closed  around 
the  z-axis  imposes  that  m be  an  integer.  We  will  also,  without  loss  of 
generality,  consider  m as  zero  or  positive  and  that  n ^ m. 
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For  h = 0,  Equation  (2. 3. 8. a)  reduces  to  the  angular  differential 
equation  obtained  in  spherical  coordinates  and  is  satisfied  by  the 


associated  Legendre  functions  of  the  first  kind  Pn(n)  and  of  the  second 


kind  Q (n)« 
n 


X-[(i  - n2)j-  P™(n)l  + [x  (0) 
dp  d r)  n mn 


m , m 


1 - n 


7]P  (n)  = 0 , (2.4.1) 

i n 


where  the  eigenvalues,  denoted  by  X (0)  when  h = 0,  are: 

mn 

X (0)  = n(n  + 1),  n > m 
mn  — 

In  that  case,  X does  not  depend  on  the  frequency  or  on  m,  yielding  a 
mn 

simple  calculation.  On  the  other  hand,  in  spheroidal  coordinates,  i.e., 

when  h ^ 0,  X depends  on  both  h and  m,  which  makes  the  computations  more 
mn 

difficult.  In  physical  problems,  we  are  interested  in  solutions  which  are 

finite  at  n = +1;  these  solutions  are  known  as  the  eigenfunctions  of  the 

first  kind  S (h,n)  of  order  m and  degree  n.  Equation  (2. 3. 8. a)  differs 
mn 

from  Equation  (2.4.1)  only  by  having  a singularity  at  infinity;  this 
suggests  an  expression  of  the  eigenfunctions  under  the  form: 


S (h.n)  = I dmn(h)Pm  (n)  . 

mn  r=o,l  r 7 m+r 


(2.4.2) 


The  prime  over  the  summation  sign  indicates  that  the  summation  is 
performed  only  on  even  values  of  r when  n-m  is  even,  and  on  odd  values  of 
r when  n-m  is  odd.  Angle  functions  of  the  second  kind  may  be  expressed 
in  an  identical  manner  from  the  associated  Legendre  functions  of  the 
second  kind,  but  their  singularities  at  n = +1  give  them  little  usefulness 
in  most  physical  problems. 

Then,  substitution  of  Equation  (2.4.2)  into  Equation  (2. 3. 8. a)  yields 
a recursion  formula  for  the  coefficients  d™n.  Use  is  then  made  of  the 
orthogonality  property  and  of  the  recursion  formula  for  the  associated 

S \ 

Legendre  functions  p (p).  The  relation  obtained  is 

n \ 


] 2 


Amdmn  + Bn>d>'in(h)  + cmdmn,(h)  = 
r r+2  r r r r-2 


(2.4.3) 


m _ (2m  + r + l)(2m  4-  r + 2)  ,2 

r = (2m  + 2r  + 3)  (2m  + 2r  + 5) 


(2. 4. 4. a) 


Dm  . , . , . . , . . 2(jn  + r)  (m  + r + 1)  - 2m  - 1 , 2 . . , . 

B = (m  + r)(m  + r + 1)  - A.  + — ~ Two — r"o — •.  t v h ,(2.4.4.b) 

r mn  (2m  + 2r  - 1)  (2m  + 2r  + 3)  ’ 


= r(r  - 1) 2 

r (2m+2r+3)(2m+2r+5) 


(2.4.4.C) 


Equation  (2.4.3)  is  a linear  homogeneous  difference  equation  of  the  second 
order.  It  must  admit  two  non-trivial  independent  solutions.  Examination 
of  the  recursion  formula  (2.4.3)  shows  that  as  r approaches  infinity 
d™n (h) /d™1^  (h)  either  increases  as  -4r^/h^  or  decreases  to  zero  as  -h^/4r^. 
Therefore,  the  latter  is  chosen  in  order  to  insure  absolute  convergence  of 


the  series  defined  in  Equation  (2.4.2).  Expressing  the  condition 

lim  dmn(h) 
r -*■  °° — = 0 

dr-2 (h) 

2 

in  Equation  (2.4.3)  leads  to  a transcendental  equation  between  X and  h . 

mn 

These  calculations  are  discussed  in  Flammer  (1957)  and  Lauchle  (1970). 

Different  methods  of  calculation  have  been  developed  to  solve  the 
very  complicated  transcendental  equation  last  obtained,  such  as  the  power 


series  developments  of  the  eigenvalues,  or  the  Bowkamp's  method  of 


approximation  (Flammer  (1975)].  The  latter  has  been  utilized  by  Van 
Buren  et  al.  (1975)  in  their  tables  of  prolate  spheroidal  angular  functions 


which  we  used  in  the  calculations  reported  later  in  this  document. 

A last  step  in  the  determination  of  the  spheroidal  wave  functions 

is  the  normalization  of  the  functions  S (h,q).  The  fact  that  S reduces 

mn  mn 

to  Pm  when  h is  zero  is  considered.  Different  methods  have  been  introduced 
n 


by  Chu  and  Stratton  (1941),  Page  (1944),  Morse  and  Feshbach  (1953), 

Meixner  and  Shafke  (1954),  and  Flammer  (1957).  For  the  same  reason  as 

In  the  preceding  paragraph,  we  will  consider  the  method  of  Meixner  and 

Shafke,  which  has  been  used  by  Van  Buren  et  al.  (1975).  The  d™n(h)  must 

be  such  that  S (h,n)  and  Pm(n)  have  the  same  normalization  factor,  i.e.: 
mn  n 

,+l  r„  /.  .+1 T m , .,2,  2 (n  + m) ! ,,  , .. 

N = / [s  (h,n)j  dn  = f [P  (n)rdn=  ~ — — r rr-  . (2.4.5) 

mn  mn  _1  n 2n  + 1 (n  - m) ! 

Substituting  Equation  (2.4.2)  into  Equation  (2.4.5),  and  using  the 
orthogonality  property  of  P™(n),  yields  the  following  normalization  for 
d”n(h): 

T (r  + ,2m)![dmn(h)]2  = ~ - /n-+-m-jr  • (2.4.6) 

r=o,l  2(m  + r)  + 1 r!  r 2n  + 1 (n  - m) ! 

Therefore,  the  coefficients  d™°  are  completely  determined,  as  well  as  the 

S (h,n). 
mn 

This  normalization  scheme  has  the  advantage  of  eliminating  the 

numerical  integration  of  the  normalization  factor  N which  is  often  used 

mn 

in  physical  applications. 

Graphs  showing  the  behavior  of  some  of  the  angular  eigenfucntions 
are  shown  in  Van  Buren  et  al.  (1975).  Different  tables  are  available  for 


values  of  d (h) , S (h,n),  A and  other  coefficients.  An  extensive 
r mn  mn 

number  of  coefficients  are  tabulated  in  Flammer  (1957)  and  eigenfunctions 
for  m = 0 and  h £ 40  are  given  in  Van  Buren  et  al.  (1975).  These  tables 
are  limited  to  low  values  of  h because  of  the  very  slow  convergence  of 
the  calculations  and  the  very  large  intermediate  values  encountered  in 
the  computations.  For  large  values  of  h,  different  approximation  formulas 
have  been  developed  by  Meixner  and  Shafke  (1954),  Silbiger  (1961),  and 


others. 
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2. 5 The  Radial  Functions 

The  radial  functions  must  satisfy  Equation  (2.3.8.b)  which  is  similar 

to  Equation  (2. 3. 8. a)  with  the  exception  of  the  parameter  range.  This 

suggests  that  we  express  the  radial  functions  R^^h.O  in  terms  of  the 

angular  functions  S (h,n)-  As  shown  in  Flammer  (1957),  solutions  of 
ran 

Equations  (2.3.8.b)  for  all  values  of  £ can  be  given  in  the  form 

R (h,£)  = /bK  (n,c)s  (h,n)dn  , (2.5.1) 

mn  a mn  mn 

where  K (n,£)  is  a Kernel  function,  discussed  below,  and  a and  b are 
m 

chosen  so  that  the  bilinear  concomittant  vanishes  at  those  limits. 

Physical  problems  usually  require  radial  functions  for  both  standing- 

wave  type  solutions  and  progressive  wave  type  solutions,  whereas  the 

angular  functions,  being  solutions  on  closed  contours,  require  only 

standing-wave  type  solutions.  So,  we  will  use  functions  of  the  first  and 

second  kind,  R^(h,£)  and  R^(h,£),  respectively,  and  functions  of  the 
mn  mn 

third  and  fourth  kind  obtained  by  complex  summation  of  the  first  two  kinds. 
As  we  have  already  noticed,  the  sphero’a'.ai  coordinate  system  degenerates 
to  the  spherical  one  as  £ becomes  infinite,  which  implies: 

, (2. 5. 2. a) 

, (2. 5.2.b) 

, (2.5.2.C) 

, (2.5.2.d) 

the  spherical  Bessel,  Neumann 


R^}(h,£)  - j (£) 
mn  _ n 

£-h® 


R^}(h,£)  -v  n (£) 
mn  _ n 

£-*» 


*<>,£)  - h<1}(£) 


and 


Rmn)(h’?)  - hn2)(C) 
mn  ^ n 

where  j (£),  n (£),  h^(£),  and  h^(£)  are 
n n n n 


and  Hankel  functions  of  the  first  and  second  kind,  respectively. 


According  to  Flammer  (1957),  the  Kernel  in  spherical  coordinates 
can  be  calculated,  and  only  the  solution  yielding  the  Bessel  functions 
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expansion  is  considered  for  our  problem,  i.e.,  after  transformation  into 

the  spheroidal  coordinate  system: 

u-  ( r\  rr\  2 m/2  2 , .m/2  ihn£  ,0  , ,v 

K (n,£)  = C(1  - q ) (5-1)  e , (2.5.3) 

m 

where  C is  an  arbitrary  constant.  Then  Equation  (2.5.1)  becomes: 

R (h.o  = c /b(l  - n2)ra/2(52  - l)ra/2eihnCs  (h,n)dn  , (2.5.4) 

mn  a mn 

and  the  bilinear  concomittant  is  given  by  [Morse  and  Feshbach  (1953)  and 

Flammer  (1957)]: 


ihnC/!  2xl+m/2r , mn 
e (1  - n ) l ( 


-o  - ihOs  (h,n)  + s'  (h,n)] 

l mn  mn 


(2.5.5) 


l - n 

Expansions  for  the  radial  functions  of  the  first,  third,  and  fourth  kind 
are  then  evaluated  by  substituting  Equation  (2.4.2)  into  Equation  (2.5.4) 
and  applying  the  three  different  sets  of  limits: 

* = -1} 

(2.5.6) 

[b  = +lj  (b  = +lj  [b  = i-J 

respectively. 

Flammer  (1957)  has  shown  that  Equation  (2.5.4)  can  be  transformed 

into: 

R (h.o  = C(52  - l)m/2r_£  , C(h)  ^2m  + ')!'  (h5)r/b’ 

mn  r-o,i  r 2m+rr ! (m  + r) ! a 


a = -1 

fa  = i°° 

(a  = -l] 

• and 

b = +1 

lb  = +lj 

b = i00 

>eihCt(i  _ t2)m+rdt 


(2.5.7) 


where  the  dmn  are  those  calculated  for  the  angle  functions.  The  evaluation 
r 

of  the  integral  in  Equation  (2.5.7)  yields  a series  of  spherical  Bessel 
and  Hankel  functions  for  the  three  sets  of  limits: 


/+1eih5ta  - t2)kdt  = 2.2^! 


^(hO 

(h£)k 


(2. 5. 8. a) 


+1  ihEt  2 k k.  ^5) 

/•+1elhU(  1 - t)kdt  = 2*k!  — : — 


( 2 . 5 . 8 . b) 


/i0°eihU(  1 - t2)kdt  = 2kk!  ^ =- 

_1  (hC) 


(2.5.8.c) 


Therefore,  Rv  '(h.O  can  be  written  as: 
mn 


R^h.O  = cii r’  d^(h)1  (2-m,  + -r-j-!  J (hO  , (2.5.9) 

mn  (hO  r=o,i  r r!  m+r 

where  the  arbitrary  constant  C is  then  determined  by  normalizing  R^^(h,£). 

mn 

A convenient  choice  is  based  on  Equation  (2. 5. 2. a)  and  the  asymptotic 

formula  for  j (£)  when  £ becomes  infinite,  i.e.: 
n 

hT  cosl^  ~ j(n  + !>J  • (2.5.10) 

mn  t,'^00  nt,  z 

After  computation  of  C,  R^(h,f,)  is  found,  according  to  the  normalization 

mn 

method  used  in  the  computation  of  the  coefficients  d™n(h)  [Equation 
(2.4.6.)]  to  be: 

■ £-;g  !-<*-;■ 


(2-m  +,  r-)-!  j ^ (h«  . 

r!  Jm+r 


(2.5.11) 


In  the  same  manner,  expressions  for  R^(h,£)  and  R^(h,£)  can  be  obtained 

mn  mn 

by  application  of  Equations  (2.5.8.b)  and  (2.5.8.c)  and  performing  the 
normalization  with  Equations  (2.5.2.c)  and  (2.5.2.d).  These  expressions 
are  given  by: 


.<2„  + r),  „(l).<»(h0  , 

r!  m+r 


(2.5.12) 
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(2) 

An  expression  for  Rmn  (h,£)  may  be  derived  using  the  classic  relation: 

e-<4)(h.«  * Rii)(h*C)  + iRin)(h*^  * (2.5.13) 

mn  mn  — inn 

as  in  Flammer  (1957),  or  by  replacing  the  spherical  Bessel  function  j 

m+r 

by  the  spherical  Neumann  functions  in  the  expression  of  R^(h,£)  [Equation 
(2.5.11)].  This  yields: 

r(2)<m)  ■ <»  - 1 1 ‘d”(h)(»  * r->'- 

mn  ’ (n  + m) ! 2 r=o,i  r v r! 

nm+r(he)  • (2.5.14) 

This  expansion  is,  in  fact,  an  asymptotic  series  which  is  not  absolutely 

convergent  for  any  finite  value  of  h£.  It  converges  only  for  large  values 

(3) 

of  h£,  and  consequently  [Equation  (2.5.13)],  the  expression  for  R (h,£) 

mn 

(4) 

and  R (h,£)  are  also  only  valid  for  large  values  of  h£,  i.e. , as 
mn 

asymptotic  forms.  Large  h£  does  not  imply  h >>  1,  £ % 1,  but  necessarily 

implies  h ^ 1,  and  £ >>  1 [Silbiger  (1961)].  An  alternate  way  of  expressing 

(2) 

R (h,£)  is  to  write  it  as  a series  of  associated  Legendre  functions  of 
mn 

both  kinds  Pm(£)  and  Qm(0,  for  £ > 1 [Hanish  et  al.  (1970)],  that  is, 
n n 

writing  that  the  radial  and  angular  functions  of  the  second  kind  are 

proportional  to  one  another  [Flammer  (1957)].  This  provides  a solution 

(2) 

for  R (h,£)  which  is  valid  for  low  values  of  h£,  and  can  be  found  in 
mn 

Flammer  (1957)  and  Hanish  et  al.  (1970). 

Tables  of  R^(h,£)  and  R^(h,£)  and  their  first  derivatives  with 
mn  mn 

respect  to  £ are  given  in  the  last  two  references,  with  a more  extensive 
range  of  the  parameters  m,  n,  h and  £ in  Hanish  et  al.  (1970). 

2.6  The  Free-Space  Green's  Function  in  Prolate  Spheroidal  Coordinates 
The  Green's  function  G(rji:’)  of  the  scalar  wave  equation  is  the 
solution  at  point  r when  a unit-strength  source  is  placed  at  point  r1 . 
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So,  G(_r|r^' ) satisfies  the  inhomogeneous  equation: 


(V2  + k2)G(r/r ' ) = 6(r  - r')  , 


(2.6.1) 


where  6 is  the  three-dimensional  Dirac  delta  function.  In  spheroidal 
coordinates  6(r_  - r ' ) is  defined  by: 


6(r  - r')  = h^h'^Vn  - n')6(C  - £')«(*  - *')  , (2.6.2) 


where  h^,  h^  and  h^  are  the  metrical  coefficients  (Flammer,  1957): 


h . 

1 - n 


(2. 6. 3. a) 


d.g2  - n2a/2 
o V.  o ’ 


£ 2 v 2 
^ 5-1 


(2.6.3.b) 


and 


ft (1  - n2)(52  - 1)]1/2 


(2.6.3.c) 


[Flammer  (1957)]. 

The  metric  coefficients  are  introduced  because  6 must  satisfy: 


/6(r  - r ' )dz ' = 1 = /2lt/+1/”[6 (r  - r ' )h'h'h’  ]dn 'd5  'd<f>  ' 

V--  O -lo,l  - - I)  U 


(2.6.4) 


with 


; ' =x+e , 


/ * a,c6(x  - x')dx'  = 1 

X =X-G 


(2.6.5) 


The  solution  of  Equation  (2.6.1)  in  the  case  when  there  are  no 


boundaries  is  called  the  free-space  Green's  function  Go(r/r')  and  is 


given  by  the  classic  formula: 


ik I r-r ' 
Go<l'l’>  " firTT^T1" 


(2.6.6) 


To  express  this  function  in  spheroidal  coordinates,  we  will  expand  it  in 
terms  of  the  natural  eigenfunctions  in  spheroidal  coordinates,  i.e.,  the 


wave  functions  previously  described.  Noting  that  G^r/r')  is  a solution 


J 


rr 
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of  the  Helmholtz  equation  for  £ / £' , and  is  symmetrical  in  £ and  £( , we 
express  it  in  the  form: 

Go(n,C,*/n\C\*’)  = m^nAmnsmn(h.n)Smn(h,n')cos  m(*  - <f>') 


R^Oj.OR^Oi,^),  £ < £' 
mn  mn 


Rli)(h»«,)R™)(h.O.  £ > £'  , 

mn  mn 


(2.6.7) 


where  represents  standing-wave  solutions,  and  progressive-wave 

mn  mn 

solutions  diverging  at  inf inity.  Therefore,  the  radiation  condition  is 

respected  for  all  observed  situations.  Furthermore,  A is  the  same  in 

mn 

both  cases  to  insure  the  continuity  of  the  Green's  function  at  £ * £'. 

We  now  integrate  Equation  (2.6.1),  applied  on  G (£/£') , over  the 
range  [£'  - e,  £'  + e] , which  yields: 

r{  (V2  + k2)G  = -/C’+V1h:1h"16(n  - n')- 

£ -e  o ^ t_e  n 

•S<£  -£')«(♦  - ♦ ')<!£  . (2.6.8) 

12  2 2 

After  multiplying  both  sides  by  (-^  d)  (£  - n ) and  making  use  of  Equation 


(2.6.5),  we  obtain: 


[<£2  - D^-  Go(n,c,(t,/n',C',4>')]J,_e 


6(n  - n,)6(<t>  — 4>’)  • (2.6.9) 


From  Equation  (2.6.7)  we  can  write: 


Jf  G(n,£,*/n’, £’,*')  = IK  •) 
of,  o m , n mn  1 


, aR^)(h,c) 

Rin>(h'e,)  » 5 < 

mn  ot. 


(2.6.10) 


with  K = A S (h,n)S  (h,n')cos  m(<J>  - $'). 
mn  mn  mn  mn 

Substituting  Equation  (2.6.10)  into  Equation  (2.6.9),  and  passing  to 
the  limit  £ = £',  the  radial  functions  and  derivatives  being  continuous 
at  £ = £'  we  get: 
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9 (n  aR™><h»5)  /vk  3R(1)(h,4) 

iU  ~ 1)«EnKmnlRm>(h*5)”J2aC R^)(h»0~- mtl  1 } - 

m,n  mn  mn  mn  9£  C=s* 

= - |«(n  - n')6(*  - *•)  . (2.6.11) 

(3) 

After  use  of  the  expression  for  (h,£)  given  in  Equation  (2.5.13),  the 
left-hand  side  of  Equation  (2.6.11)  becomes: 

, m aR^^h.O  m 3R(2)  (h,£) 

ur  - lUKtR^Oi.c)-^ + iRin  (fc.O-55^ 

m,n  mn  mn  3£,  mn  3£ 

aR^)(h,o  m SR^^h.c) 

- «^)<».e)-asR lR»n  — •!»«•  • 

Then,  using  the  Wronskian  relation  [Flammer  (1957)]: 

(n  3Rfn)(h’C)  (o\  3Rli)(h»^)  I 

Rm „ (h,C) tz R (n,£) rz = ^ 

mn  3£  mn  3£ 


h(r  - 1) 

in  this  last  expression  yields  a new  form  for  Equation  (2.6.11),  i.e.  : 

~ EK  = - |«(n  - n ’ ) <5  (<J>  - *’) 
h m , n mn  d 


or 

T.  A S (h,n)s  (h,n’)cos  m(4>  - $’)  = ik6(n  - n')6(4>  - $')  . (2.6.12) 
m,n  mn  mn  mn 

We  now  multiply  both  sides  of  Equation  (2.6.12)  by  Sp^(h,n)cos  p<)>,  and 
integrate  over  the  range  of  n and  <J>: 

/+1/2its  (h,q)cos  p^  Z A S (h,n)S  (h,n')cos  m(«]i  - <J> ' ) d d>d n * 

”1  o m,n  mn  mn  mn 

- /+1/2l,S  (h,n)cos  p<(>  ik5(n  - n')6(*  - ♦,)d*dn  • (2.6.13) 

“1  o P9 
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where  N is  the  norm  of  S (h,n')  defined  in  Equation  (2.4.5), 

pq  pq 

is  the  Kronecker  symbol. 


Therefore,  we  have  determined  the  value  of  the  constant  A 

pq 

ike 

pq  2„Npq 


and  6 

op 


(2.6.15) 


with 


e = 2 - 6 . (2.6.16) 

P op 

The  arbitrary  choice  of  p and  q in  Equation  (2.6.13)  makes  it  possible  to 
apply  Equation  (2.6.15)  for  all  possible  values  of  the  set  (p,q). 

Substituting  Equation  (2.6.15)  into  Equation  (2.6.7)  yields  the 
final  expression  for  the  free  space  Green's  function: 


Go(n,C,4/n\£\<f ) 


ik  " 

E e 

2tt  m»o  m n=m 


s (h.n)S  (h,n' ) 

ran  mn 


mn 


cos  m(4>  — <t>  * ) 


R^)(h,C)R^)(h,e')  , C < C' 

mn  mn 

R^^h.c'jR^^h.c)  , e > v 

mn  mn 


(2.6.17) 
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CHAPTER  III 


ANALYSIS 


3. 1 Introduction 

In  Chapter  II,  we  developed  expressions  for  the  spheroidal  wave 
functions  and  free-space  Green's  function.  They  are  the  tools  needed 
to  solve  the  problems  that  will  be  presented  here.  Our  purpose  is  to 
calculate  the  ratio  of  the  scattered  pressure  to  the  directly  incident 
pressure  far  away  from  a prolate  spheroid  which  is  insonified  at  a given 
frequency  by  a point  source  located  at  a finite  distance  on  the  axis  of 
the  spheroid  (Figure  2).  Two  types  of  boundary  conditions  on  the  spheroid 
are  considered:  a perfectly  rigid  surface,  and  a perfect  pressure  release 
surface.  Furthermore,  combinations  of  two  sources  of  the  same  strength 
are  investigated,  i.e.,  two  sources  in  phase  and  two  incoherent  sources. 

In  this  section,  we  derive  the  expressions  for  the  reflection  factor 
corresponding  to  those  different  situations. 


3.2  Case  of  a Prolate  Spheroid  with  Perfectly  Rigid  Surface 
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As  we  will  treat  the  problem  of  infinite  r,  we  have  £ > Using  the 

free  space  Green's  function  given  by  Equation  (2.6.17)  for  this  limit,  we 
can  express  the  incident  pressure  by 

„ S (h,n)S  (h.l) 

— — GO  n ' 


P = p Z - 
t o n=o 


on 


on 


N 


on 


(3.2.1) 


with 


k pcQ 


(3.2.2) 


We  now  want  to  express  the  pressure  at  the  observation  point,  which  is  due 
to  the  scattering  of  sound  by  the  spheroid  when  insonified  by  the  point 
source.  We  call  this  pressure  the  scattered  pressure  and  denote  it  by  p^. 
It  is  expedient  to  express  pgl  in  terms  of  the  spheroidal  wave  functions, 
i.e. : 

Rsl  ■ n£oVon(h’'1>Ron)(h-£>  • «-2-3> 

The  radial  function  of  the  third  kind  is  chosen  because  it  represents  a 
progressive  wave  diverging  at  infinity,  which  is  the  required  behavior. 

The  boundary  condition  for  the  rigid-type  surface  is  given  by  the 
Neumann  boundary  condition: 


3p.  8p  , 


= 0 at  ( ■ ( and  £ 
o o 


(3, 2. A) 


35  35 

where  p,  and  p . are  the  incident  and  scattered  pressures,  respectively, 
1,0  si , o 

on  the  boundary  surface.  On  the  surface,  the  incident  pressure  is  written 
as  [Equation  2.6.17)]: 

- s (h,n)s  (h,n’> 


p.  = p £ 

*1,0  ro  n»o 


on 


on 


on 


R^)(h,C)R^)(h,C’)  at  e = 5 , (3.2.5) 

on  on  o 


and  Equation  (3.2.4)  can  be  written  as: 
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P E 
o n=o 


S (h,n)S  (h,n')  / -i  v 

R2  ,(h>5„)E  (h^'J  B " E AS  (h,n) 

O n q Of!  nsn  n nn  * 


N 


on 


n=o  n on 


R£),(h’Sj  * (3.2.6) 
on  o 


where 


and 


RL  <h-V  ■ u k-( 


n)  3R<3)(h,s) 

Ro„  <h-V  - U-E 


(3. 2. 7. a) 


(3 . 2 . 7 .b) 


The  boundary  condition  must  be  respected  on  the  whole  surface  of  the  body 
and  for  any  source  location  which  means  that  Equation  (3.2.6)  must  be  true 
for  all  values  of  n and  C'  (provided  £'  > CQ) • Therefore,  Equation  (3.2.6) 
must  be  verified  separately  for  each  order  n.  Hence,  after  division  by 


-S  (h,n)  we  obtain: 
on 


-P  °n*  ' ~ R(1),(h  ,r,  )R(3)(h,£')  = A R(3) ' (h,C  ) 
o N on  o on  n on  o 

on 


(3.2.8) 


for  all  values  of  n.  This  gives  the  expression  for  the  coefficients  A^. 
After  substituting  this  value  into  Equation  (3.2.3),  we  obtain  the 
expression  of  Pg^: 

S (h,n)S  (h,l)  R(1),(h,£  ) ...  ... 

— “ on  ’on  on  o „(3)„  .,,„(3)„ 

Psl  = -po  n=o  N- ToyrrTT  Ron  (h’5  )Ron  (h’?)  * (3-2’9) 


on 


R™,(h’S0> 


From  the  expression  for  p^  and  p , we  can  now  write  the  reflection  factor 


R1(h,C,n)  = . 


(3.2.10) 


As  we  are  interested  in  the  far-field  reflection  factor,  we  will  replace 

.O), 


(h,£)  by  its  asymptotic  expansion: 
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R(3)(h,£)  - i eih5e-i  2 (n+1)  = (-i)n+1 

on  ^ h£  H 


(3.2.11) 


so  that  the  final  expression  for  the  ref lection  factor  at  Infinity  in  the 
case  of  the  rigid  spheroid  is  given  by: 


R(1)'(h,£  ) ... 

r a°/l  s on  o „(3)„  r(w  ,sn 

Z A (h.n)-^ R (h,£')(-i) 

no  n r,  ( j)  i /i  v on 


R^Oi.o)  = -- 


Rnn',(h**„> 
on  o 


f A*(h,n)R(1)(h,£')(-i)n 

n=o  n on 


(3.2.12) 


with 


s (h,£)s  (h,i) 

iO/i  \ on  on 

A (h,n)  = n 

n N 


(3.2.13) 


on 


We  notice  that  this  result  is  valid  for  any  value  of  the  source  angular 
coordinate  n'  provided  that  symmetry  with  respect  to  the  z-axis  is 
conserved.  A value  of  n'  different  from  +1  or  -1  would  then  describe  a 
problem  with  a ring  source  defined  by  the  intersection  of  the  spheorid  £' 
and  the  hyperboloid  defined  by  n' . 


3 . 3 Case  of  a Prolate  Spheroid  with  Pressure  Release  Type  Surface 

The  situation  is  analogous  to  the  previous  one  except  for  the  type 
of  boundary  condition.  The  expression  for  the  incident  pressure  remains 
unchanged,  and  we  still  express  the  scattered  pressure  (denoted  by  p^) 
in  the  form  given  in  Equation  (3.2.3). 

Then,  we  write  a new  boundary  condition,  which  is  a Dirichlet  boundary 
condition: 

(3.3.1) 


p + p . = 0 at  C = ( and  £ < £ ' 

*i,o  s2 ,o  o o 


Keeping  the  same  notation  as  in  the  previous  section,  we  find: 


1 

J 
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Snn<h^>S  (h.D 


p z 

* o n=o  N 


Da>,u , x»o) 


- K'  (h,Or'(h,?')  = -I  A S (h,n) 

on  o on  n=o  n on 


Ro»)(h^o> 


(3.3.2) 


which  determines  the  value  of  for  all  n.  Hence,  the  expression  for  the 
scattered  pressure  p ^ is 

Js2  - ‘Pc  Ro»)<h'{')Eon><h-t>  ' <3-3-3> 

on  o 

Making  the  same  approximation  as  in  Equation  (3.2.11)  we  obtain  the 
expression  for  the  reflection  factor  at  infinity  for  the  case  of  the 


spehroid  with  pressure  release  surface: 


R2(h,r,) 


„ R(1)(h,£  ) ... 

nVOh’n)~?T) ~ R n )(h,C')(-i)n 


(3.3.4) 


3.4  Combination  of  Two  Sources 

We  now  consider  the  combination  of  two  sources  of  equal  strength  on 
the  axis,  at  finite  distances  from  the  spheroid.  Two  cases  are  investigated: 
the  sources  are  in  phase,  and  the  sources  are  incoherent. 


In  the  first  case,  we  simply  apply  the  superposition  theorem  by  adding 
the  incident  pressures  due  to  each  source  on  one  hand,  and  the  scattered 
pressures  due  to  each  source  on  the  other  hand.  This  can  be  written  as: 


Pi  = Pi,l  + Pi , 2 


(3. 4.1. a) 


Ps  = ps,l  + Ps , 2 * 


(3.4.1.b) 


where  subscript  1 and  2 refer  to  Sources  1 and  2,  respectively,  and  where 


Pj  is  the  resulting  incident  pressure,  and  p the  resulting  scattered 
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pressure.  We  will  denote  scattered  pressures  p . , p and  p in  the 

Si  Si,  I S.L,  £ 

rigid  type  boundary  problem,  and  pg2,  pg2  ^ and  pg2  2 in  the  pressure 
release  type  boundary  problem.  Generally,  we  can  write  the  reflection 
factor  in  the  form: 


ps  Ps  i + ps  2 

R(h,n)  -zrm  _ — 

Pi  pi,l+pi,2 


(3.4.2) 


and  the  computation  will  be  made  accordingly. 

For  the  incoherent  sources,  statistical  theory  shows  that  the  energies 
add  and  not  the  pressures.  In  terms  of  pressure,  this  can  be  expressed  by 
adding  the  square  of  the  magnitude  of  the  pressures.  So  we  obtain: 

i — i 2 i—  i 2 . i — i 2 


lpi,l'  + ' pi , 2 1 


(3. 4. 3. a) 
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(3.2.2).  Therefore,  we  can  apply  the  reciprocity  property  of  the  Green's 
functions  to  our  problem,  in  the  same  manner  as  Horse  (1968)  applied  it 
to  the  problem  of  diffraction  by  a cylinder.  We  can  say  that  the  incident 
pressure  p^  at  the  point  of  coordinates  (+°°,n),  due  to  a point  source 
located  at  (£',n'),  is  the  same  as  the  pressure  created  at  (£*,n')  by  the 
same  source  located  at  (+<®,ri).  The  same  reciprocity  theorem  can  be  applied 
to  the  scattered  pressure.  Therefore,  it  is  also  true  for  the  reflection 
factor. 

As  we  have  considered  £ the  wave  impinging  on  the  body  in  the 

reciprocal  situation  is,  in  fact,  a plane  wave  which  is  incident  under 
the  aspect  angle  0 = cos  ^ n.  Consequently,  our  reflection  factor  is  also 
the  reflection  factor  at  the  observation  point  on  the  z axis  defined  by 
( £ ' , +1 ) when  the  body  is  insonified  by  a plane  wave  impinging  under  the 
aspect  angle  0 = cos  ^ n with  respect  to  the  positive  z semi-axis. 


CHAPTER  IV 


NUMERICAL  COMPUTATIONS  AND  RESULTS 


4 . 1 Introduction 

In  this  section,  we  present  the  calculations  which  have  been  carried 
out  according  to  the  formulas  derived  in  the  last  chapter . The  computations 
have  been  made  on  an  IBM  370/167  digital  computer  located  at  The 
Pennsylvania  State  University  campus,  using  the  WATFIV  compiler  which 
allows  computations  directly  in  complex  variables  as  expressed  in 
Equations  (3.2.12),  (3.3.4),  (3.4.2),  and  (3.4.4),  Single  precision 
was  used  in  all  of  the  computations. 

4 . 2 Range  of  the  Parameters 

As  seen  in  the  previous  chapter,  we  have  considered  two  types  of 

boundary  conditions  and  three  source  situations.  The  body  radial  coordinate 

c,  remained  constant  in  all  the  calculations,  i.e.: 
o 

5 = 1.01 

o 

which  corresponds  to  a fineness  ratio  [Equation  (2.2.4)]: 


L 


In  the  case  of  a single  source,  ten  different  locations  have  been  considered 
along  the  positive  z axis:  1.02  £ £'  _<  1.20  in  increments  of  0.02. 

The  values  of  the  radial  wave  functions  of  the  first  and  second  kind 
and  their  first  derivatives  have  been  taken  from  the  tables  by  Hanish  et 
al.  (1970).  Then  the  function  of  the  third  kind  was  written  according 


to  Equation  (2.5.13). 
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In  the  case  of  two  sources,  we  considered  one  source  fixed  at  the 
position  £'  = 1.02,  the  second  source  successively  taking  the  nine  remaining 
positions  denoted  in  the  range  above. 

For  all  source  locations,  we  determined  the  directivity  pattern  of 
the  reflection  factor  for  three  values  of  the  parameter  h: 

h = 20,  35,  40  . 

Those  calculations  were  made  for  values  of  the  aspect  angle  0 = cos 
from  0 to  180  degrees,  in  increments  of  two  degrees.  The  values  of  the 
angular  wave  functions  for  0 £ 90°  were  obtained  from  the  tables  of  Van 
Buren  et  al.  (1975).  For  0 > 90°,  we  used  the  symmetry  relation: 

S (h,-n)  = (-l)nS  (h,n)  (4.2.1) 

on  on 

which  can  also  be  written: 

S (h,n  - 0)  = (-l)nS  (h,0)  . (4.2.2) 

on  on 

For  all  of  the  source  locations,  we  also  determined  the  frequency 
response  of  the  reflection  factor  for  three  values  of  the  aspect  angle: 

0=0,  90,  180°  . 

Those  calculations  were  carried  out  for  the  following  thirteen  values  of 
h: 

h = 0.1,  0.4,  0.6,  1.0,  2.0,  4.0,  6.0,  10.0,  12.0,  14.0,  16.0,  18.0,  20.0. 

We  have  listed  in  Tables  I,  II,  and  III  the  various  parameters 
corresponding  to  the  different  source  locations  and  the  nondimen sional 
frequencies  h = 20,  35,  and  40.  In  those  tables,  d is  the  interfocal 
distance,  Lq  the  length,  and  Rq  the  minor  radius  of  the  spheroid,  X the 
wavelength,  i'  the  distance  from  the  single  source  located  at  to  the 
closest  pole  of  the  spheroid,  £"  the  radial  coordinate  of  the  second 
source  and  t,"  its  distance  from  the  pole  of  the  spheroid,  and  a the 
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separation  distance  between  the  two  sources.  Those  values  were  determined 
from  the  following  simple  relations: 


v = -a1  - c ) , 

2 O 


2n  d 

h = X 2 ’ 


(4. 2. 3. a) 
(4.2.3.b) 


R = y42  - 1 , 

o 2 o 


a = |(C"  - £')  = A"  - *' 


(4.2.3.C) 

(4.2.3.d) 


and 

L * 5 d . (4.2.3.e) 

o o 

We  note  that,  because  is  very  close  to  1,  Lq  = £Qd  is  not  very 
different  from  d. 

As  the  tables  give  values  of  the  angular  and  the  radial  functions 

only  for  n covering  the  range  from  0 to  49,  an  important  point  was  to 

check  the  convergence  of  the  three  series  used  for  the  computation  of 

p./p  f p , /p  , and  p „/p  . We  have  checked  the  convergence  of  these 
ri  o si  o s2  o 

series  as  a function  of  the  upper  limit  of  summation  N,  for  h = 20,  35, 
and  40.  Some  of  the  results  are  presented  in  Figure  3.  In  general,  it 
was  found  that  a number  of  terms  between  h/2  and  h is  necessary  for 
convergence,  in  the  range  £'  = 1.02  to  1.20.  Therefore,  we  used  50  terms 
in  the  computations  of  directivity  patterns  at  h = 20,  35,  and  40,  and 
25  terms  for  frequency  response  curves  between  h = 0 and  h = 20. 


4.3  Results 

In  spite  of  the  fact  that  all  the  calculations  have  been  made  for 
both  the  rigid-type  and  the  pressure  release-type  boundaries,  the  bulk 


of  the  results  presented  are  for  the  former  case;  only  a few  typical 
results  are  presented  for  the  latter  case. 
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4.3.1  Case  of  the  rigid  boundary.  For  the  rigid  spheroid  insonified 
by  both  one  and  two  sources,  we  have  plotted  directivity  patterns  showing 
the  modulus  of  the  reflection  factor,  j R | as  a function  of  9 in  Figure  4 
through  20  and  the  variations  of  |r|  with  h are  shown  in  Figure  21  through 
35.  In  all  directivity  curves,  the  source(s)  is  (are)  located  on  the  0 = 
0°  axis. 


In  the  case  of  a single  source,  the  directivity  patterns  (Figures  4- 
9)  are  very  similar  to  each  other.  They  exhibit  a nearly  omnidirectional 
pattern  for  0 _<  90°,  while  for  0 > 90°,  a lobe  develops  in  the  forward 
direction.  For  a given  value  of  h,  successive  curves  are  similar  to  each 
other  when  £'  is  larger.  In  particular,  curves  obtained  for  * 1.02  are 
only  slightly  different  from  the  curves  for  > 1.02.  When  the  separation 
distance  l'  increases,  the  overall  level  of  the  reflection  factor  is  seen 
to  decrease.  This  may  be  due  to  the  lower  amount  of  energy  impinging  on 
the  spheroid  (spherical  spreading  loss).  But  with  £ * / A remaining  constant, 
the  level  of  [ R | is  higher  when  a higher  frequency  in  considered,  which 
seems  to  imply  diffraction  phenomena  becomes  more  important. 

In  the  case  of  the  combination  of  two  in-phase  sources  (Figures  10- 
14),  the  directivity  patterns  differ  from  the  patterns  obtained  with  a 
single  source,  and  are  more  irregular.  Of  particular  interest  is  the 
formation  of  cusps  at  some  angles.  These  cusps  should  not  be  interpreted 
as  a sudden  increase  in  scattered  pressure  as  this  is  not  the  case.  They 
result  because  of  a null  in  the  incident  pressure  field  defined  by  Equation 
(3. 4. 1.2).  The  two  in-phase  sources  separated  by  a distance  a [see  Table 
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II  and  Equation  (4.2.3.d)]  represent  a doublet  whose  directivity  function 
is  given  by 

D,  = cos(~  cos  0)  = cos(rp  cos  6)  . (4.3.1) 

doublet  2 A 

The  nulls  in  the  incident  pressure  field  (which  correspond  to  cusps  in 
the  reflection  factor),  therefore,  occur  at  angles  for  which 

f COS0C  = ± » (4.3.2) 

where  p is  a positive  integer.  It  further  follows  that  cusps  will  not 

occur  if  a/A  < y. 

The  values  of  a/A  for  the  various  source  locations  and  reduced 
frequencies  are  tabulated  in  Table  III.  Also  shown  in  this  table  are  the 
critical  angles  (6^)  calculated  from  Equation  (4.3.2);  obviously,  angles 
(tt-0c)  are  also  critical;  they  correspond  to  the  minus  sign  in  the  right- 
hand  side  of  Equation  (4.3.2).  Clearly,  these  angles  correspond  closely 
to  the  aspect  angles  of  reflection  factor  where  a cusp  occurs  (Figures  11-14). 

In  the  case  of  two  closely-spaced  incoherent  source  (Figures  16-20), 
the  directivity  pattern  for  the  incident  pressure  in  the  far  field  does 
not  depend  any  longer  on  the  phase  relationship  between  the  source,  i.e., 
|pil  is  essentially  omnidirectional.  For  the  scattered  pressure,  there 
is  a simple  addition  of  energies  and  therefore  an  "averaging"  of  the  two 
scattered  fields.  Because  of  the  spherical  spreading  of  the  wave  incident 
on  the  body,  the  sound  wave  due  to  the  outermost  source  is  weaker  than 
that  due  to  the  innermost  source.  As  a consequence,  the  resulting  ref  lection 
factor  is  very  close  to  that  one  due  to  the  single  source  located  at  = 
1.02,  for  a given  frequency,  i.e.,  | R ( changes  little  as  £"  increases. 

We  could  expect  that  for  larger  values  of  £",  the  energy  due  to  the  outermost 
source  would  be  negligible,  and  the  resulting  reflection  factor  would  be 
identical  to  the  reflection  factor  due  to  the  innermost  source  alone. 
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The  frequency  response  curves  of  |r|  (Figures  21-35)  also  show  the 
similarities  and  difference  cited  above. 

In  the  case  of  a simple  source  (Figures  21-25),  the  level  of  |r| 
decreases  over  the  entire  frequency  range  (0.1  < h < 20)  as  £'  increases. 

In  the  case  of  two  sources  in  phase  (Figures  26-30),  especially  for 
£"  > 1.10,  the  frequency  response  curves  are  seen  to  be  less  regular. 

This  is  due  to  the  complexity  of  the  incident  field. 

In  the  case  of  two  incoherent  sources  (Figures  31-35),  all  frequency 
response  curves  are  similar  to  the  frequency  response  curve  of  a single 
source  at  £'  = 1.02. 

A generality  for  all  of  the  curves  is  that  at  low  frequencies,  i.e. 
h < 1.5,  | R | increases  with  a slope  of  6 dB/octave  in  the  9=0°  direction 
and  12  dB/octave  in  the  0 = 90°  direction.  In  the  direction  9 = 180°, 

| R | increases  linearly  with  a slope  slightly  less  than  6 dB/octave  over 
the  frequency  range  of  computation.  We  note  that  for  h < 1.5,  all  the 
curves  have  the  same  behavior  in  a given  direction  (0=0°,  9 = 90°  or 
0 = 180°),  but  also  that  for  h < 1.5  and  £'  not  greater  than  1.20  we 
have: 

X > 2d 

and 

V < 0. 2d  . 

This  means  that  all  the  points  on  the  surface  of  the  spheroid  are  almost 
in  phase  and  the  source  is  relatively  close  to  its  pole;  therefore,  the 
situation  is  similar  to  the  case  when  the  source  is  on  the  body  itself. 

It  is  noted  that  the  frequency  response  curves  begin  to  deviate  from  a 
constant-slope  character  for  reduced  frequencies  greater  than  approximately 
1.5.  This  is  to  be  expected  because  d > X/2  for  this  range. 


Using  Table  IV,  we  have  expressed  |r|  as  a function  of  the  nondimens ional 
parameter  ki',  in  the  case  of  a single  source  and  for:  0=0°,  90°,  180°, 
and  h = 16,  20,  35,  40.  Figure  36  shows  these  results.  It  appears  that 
kJl'  is  not  a perfect  scaling  parameter.  However,  based  on  these  curves,  it 
does  appear  that  the  magnitude  of  the  reflection  factor  uniformly  decreases 
as  the  source  moves  farther  away  from  the  pole  of  the  spheroid. 

4.3.2  Case  of  the  pressure  release  boundary.  As  mentioned  before, 
we  present  only  a few  of  the  results  for  a pressure  release  boundary. 

Curves  are  presented  for  two  source  locations:  £’  = 1,02,  £'  «=  1,12,  and 
their  combinations  (Figures  37  to  41). 

The  directivity  patterns  for  a single  source  are  very  regular,  having 
a slightly  higher  level  in  the  forward-scattered  direction.  This  regular ity 
is  due  to  the  type  of  boundary  condition.  Indeed,  the  source  is  always 
relatively  close  to  the  body  where  we  know  that  for  the  source  on  the  body 
itself,  an  omnidirectional  pattern  always  results  with  such  a houndary 
condition.  In  comparing  the  calculations  of  the  soft  spheroid  with  those 
of  the  rigid  boundary,  for  a common  set  of  parameters,  the  overall  level 
of  the  reflection  factor  is  seen  to  be  slightly  less,  indicating  that 
diffraction  is  not  as  important. 

The  frequency  response  curves  show  that  |R|  remains  almost  constant 
for  low  frequencies  (h  < 1,0),  This  is  to  be  expected  hecause  of  the  type 
of  boundary  condition.  Because  of  the  spherical  spreading  of  the  impinging 
wave  as  h or  £'  increases,  the  level  of  |r|  decreases  and  the  situation 
becomes  different  from  the  limiting  case  of  the  source  on  the  hody.  This 
results  in  a shorter  range  of  h for  which  the  low-frequency  behavior  is 


linear . 
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CHAPTER  V 

SUMMARY  AND  CONCLUSIONS 

A theoretical  investigation  of  the  scattering  of  sound  by  a prolate 
spheroid  insonified  by  a point  source  or  a combination  of  point  sources  on 
the  axis  has  been  performed.  Two  limiting  types  of  boundary  conditions 
on  the  body  have  been  considered,  i.e. , Neumann  and  Dirichlet  boundary 
conditions.  Both  directivity  patterns  and  frequency  response  curves  have 
been  obtained.  This  work  has  been  limited  to  low  frequencies  (h  _<  40) 
where  numerical  computations  have  been  performed  using  tables  of  the 
spheroidal  wave  functions. 

For  a single  source  located  on  the  axis  of  symmetry  at  various  positions 
relative  to  the  pole  of  the  rigid  spheroid,  the  directivity  characteristics 
were  found  to  be  essentially  omnidirectional.  When  two  in-phase  sources 
were  placed  on  the  axis,  the  directivity  characteristics  were  found  to 
contain  much  variation  with  aspect  angle.  If  the  two  sources  were  assumed 
incoherently  related,  the  patterns  approached  the  shape  of  the  corresponding 
pattern  for  the  single  source  closest  to  the  pole,  particularly  when  the 
second  source  was  relatively  far  from  the  first.  The  reflection  factor 
was  found  to  increase  with  frequency  by  6 dB/octave  for  the  forward  and 
backscattered  directions,  and  by  12  dB/octave  for  the  broadside  aspect 
angle. 

As  this  work  has  been  strictly  theoretical , it  would  be  of  great  interest 
to  obtain  experimental  results,  particularly  for  the  source  positions  close 
to  the  body,  and  for  the  combinations  of  incoherent  sources.  Emphasis 
could  be  placed  on  frequency  responses  which  have  not  been  commonly 
investigated . 

ji 
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It  also  seems  appropriate  to  extend  the  present  results  to  higher 
frequencies.  As  tabulated  values  of  the  eigenfunctions  are  not  available 
for  h > 40,  it  would  be  necessary  to  use  asymptotic  expansions  as  given 
in  Silbiger  (1961)  and  Slepian  (1965),  for  example.  Those  theoretical 
results  could  also  be  compared  to  experimental  results  obtained  at  higher 


frequencies . 
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Table  I 

51  * /d  and  fc'/R  for  Different  Source  Locations 
o 


^0 

C 

i'/d 

l'/R 

o 

1.01 

1.02 

0.005 

0.035 

- 

1.04 

0.015 

0.106 

- 

1.06 

0.025 

0.176 

- 

1.08 

0.035 

0.247 

- 

1.10 

0.045 

0.317 

- 

1.12 

0.055 

0.388 

- 

1.14 

0.065 

0.458 

- 

1.16 

0.075 

0.529 

- 

1.18 

0.085 

0.600 

— 

1.20 

0.095 

0.670 
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Table  II 

a/d  and  a/R  for  Different  Sets  of  Source  Locations 
o 


e" 

a/d 

a/R 

o 

1.02 

1.04 

0.010 

0.071 

- 

1.06 

0.020 

0.141 

- 

1.08 

0.030 

0.212 

- 

1.10 

0.040 

0.282 

- 

1.12 

0.050 

0.353 

- 

1.14 

0.060 

0.423 

- 

1.16 

0.070 

0.494 

- 

1.18 

0.080 

0.564 

_ 

1.20 

0.090 

0.634 
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Figure  3.  Convergence  of  the  Incident:  and  Scattered  Pressure  Expansions; 
h = 40;  V = 1.12 
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Figure  9.  Directivity  Patterns,  Rigid  Case,  One  Source;  h 
1.10,  1.14,  1.16,  ard  1.20 


Figure  10.  Directivity  Patterns,  Rigid  Case,  Two  Sources  in  Phase; 

1.02;  h = 20;  C = 1.04,  1.06,  1.08,  1.10,  1.12,  and  1.14 
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Figure  11.  Directivity  Patterns,  Rigid  Case,  Two  Sources  in  Phase;  £" 
1.02;  h = 20;  £"  - 1.16  and  1.20;  h - 35;  * 1-04,  1.06 

1.08,  and  1.10 


Figure  12.  Directivity  Patterns,  Rigid  Case,  Two  Sources  in  Phase; 

1.02;  h = 35;  C = 1.14;  h = 40;  ?"  = 1.04,  1.06,  1.08  and 


Figure  16.  Directivity  Patterns,  Rigid  Case,  Tvo  Incoherent  Sources;  t 
1.02;  h * 20;  C m 1.08,  1.12,  1.16,  and  1.20 


Figure  17.  Directivity  Patterns,  Rigid  Case,  Two  Incoherent  Sources;  £ 
1.02;  h «=  35;  £"  = 1.04,  1.06,  1.10,  1.14,  and  1.18 
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Figure  21,  Frequency  Response  Curves,  Rigid  Case,  One  Source;  (a)  £ 
1.02,  (b)  = 1,04 
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Pigure  23.  Frequency  Response  Curves,  Rigid  Case,  One  Source;  (a)  £ 
1.10;  (b)  5»  = 1.12 
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Figure  36.  Reflection  Factor  as  a Function  of  ki.',  Rigid  Case,  One 
Source;  (a)  9=0*;  (b)  0 = 90°;  (c)  0 = 180 
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